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ABSTRACT  OF  WORK  UNDER  ARO  CONTRACT  DAAL03-86-K-0147 

The  project,  "Spectroscopic  Investigations  of  nonambient  systems  with  Millimeter  and 
Submillimeter  Probes,"  addressed  a  broad  range  of  scientific  topics  and  technological 
developments  important  to  this  spectral  region.  For  the  purposes  of  this  discussion,  we  have 
divided  the  report  on  research  results  into  three  parts:  (1)  Energy  transfer  in  polyatomic 
molecules,  especially  those  of  importance  in  far  infrared  lasers;  (2)  Molecular  collisional 
processes  at  very  low  temperature  (<  5  K),  and  (3)  Development  and  uses  of  millimeter  and 
submillimeter  technology.  Also  included  is  a  list  of  personnel  who  worked  on  these 
projects. 

I.  INTRODUCTION  AND  BACKGROUND 

This  project  addressed  the  development  and  use  of  spectroscopic  techniques  for  the 
study  of  molecular  systems  under  nonambient  conditions.  Emphasis  was  placed  on  the  use 
of  millimeter  and  submillimeter  spectroscopic  methods  as  diagnostic  probes  of  these  environ¬ 
ments  although  many  of  these  techniques  are  also  applicable  in  other  spectral  regions.  In 
some  cases,  the  mm/submm  spectral  region  was  especially  advantageous  for  the  work  while 
in  others  it  was  merely  convenient  because  of  our  experience  and  equipment  base  in  this  part 
of  the  spectrum. 

Much  of  what  is  proposed  is  based  upon  previous  developments  in  our  laboratory.  In 
this  Introduction,  we  will  briefly  discuss  some  of  the  basis  for  this  work;  in  later  sections 
more  of  the  details  and  references  are  provided. 

A.  The  mm/submm  Spectral  Region 

1.  The  spectroscopic  technique:  Ti  e  mm/submm  spectroscopic  technique  which  was 
used  in  many  of  these  experiments  has  been  developed  over  a  period  of  years.1'3  Briefly,  the 
nonlinear  harmonic  generators  which  produce  the  mm/submm  radiation  are  driven  by  klys¬ 
trons  or  traveling  wave  tube  amplifiers4  in  the  region  around  50  GHz.  These  microwave 
sources  are  in  turn  referenced  to  or  driven  by  sources  which  are  controlled  by  frequency  syn- 

'  ■  k 

thesizers,  effectively  giving  a  synthesized  source  which  is  conveniently  and  continuously 
tunable  throughout  the  region  between  ~  100  •  1000  GHz.  This  mm/submm  energy  is 
radiated  quasi-optically  through  the  atomic  and  molecular  system  being  studied  and  detected 
by  an  InSb  detector  operating  at  1.5  K  or  a  0.3  K  germanium  bolometer.  Associated 
electronics  are  used  for  frequency  measurement  and  signal  recovery.  Systems  of  this  type 
have  proven  to  be  reliable,  easy  to  operate,  and  reasonably  inexpensive.  Because  of  these 
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attributes,  a  number  of  similar  systems  have  now  been  built  at  laboratories  around  the  world 
and  have  become  the  standard  for  a  wide  range  of  molecular  studies  in  this  spectral  region. 

2.  Characteristics  of  mm/submm  spectroscopy:  Much  of  the  work  discussed  in  this 
report  depended  upon  several  important  attributes  of  our  technique  for  spectroscopic  studies 
in  the  mm/submm  spectral  region.  Among  these  are: 

a)  It  is  a  very  sensitive  technique.  Under  equilibrium  conditions,  rotational  absorp¬ 
tion  coefficients  increase  as  v2-»  v3  and  peak  typically  in  the  vicinity  of  1000  GHz.  Under 
optimum  non-equilibrium  conditions,  achievable  in  our  pumped  experiments,  concentrations 
as  small  as  100  cm*3  (~10'14  Torr  partial  pressure)  can  be  studied  with  modest  signal  averag¬ 
ing.  Calculations  based  on  a  more  aggressive  set  of  assumptions  yield  concentrations  as 
small  as  1  cm'3. 

b)  Because  of  the  broad  tunable  coverage  in  the  frequency  regime  that  corresponds 
to  many  of  the  rotational  transitions  of  the  small,  fundamental  species  that  we  wish  to  study, 
transitions  can  be  chosen  on  the  basis  of  their  diagnostic  value  rather  than  to  satisfy  coinci¬ 
dence  criteria.  Furthermore,  because  of  the  high  resolution,  spectral  overlap  is  very  rare. 

c)  It  is  an  absorption  technique.  Because  the  transitions  that  are  observed  are  rota¬ 
tional  transitions  with  accurately  known  transition  moments,  the  absolute  absorption  data  can 
be  accurately  converted  into  molecular  information. 

d)  Because  we  directly  observe  rotational  transitions,  the  method  is  especially 
sensitive  to  rotational  nonequilibria.  For  example,  at  300  K  and  300  GHz  a  5  %  change  in 
the  population  of  one  of  the  states  involved  in  a  transition  produces  a  100  %  change  in  the 
observed  line  strength. 

e)  Because  this  is  the  frequency  regime  in  which  a  very  large  number  of  optically 
pumped  lasers  have  been  discovered,  we  can  take  advantage  of  the  hundreds  of  man  years 
that  have  been  spent  searching  for  C02  laser  pump  coincidences. 
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B.  Applications 

Over  a  number  of  years  we  have  used  these  mm/subram  techniques  for  a  wide  variety  of 
spectroscopic  investigations.  In  the  sections  below  some  of  those  that  are  related  to  the 
reported  work  are  briefly  discussed. 

1.  Spectroscopy  of  small,  fundamental  molecular  species:  The  principal  initial  applica¬ 
tion  of  this  spectroscopic  technique  was  its  use  to  develop  a  basic  spectroscopic  understand¬ 
ing  of  small,  fundamental  molecules  in  the  mm/submm.  Particular  emphasis  was  placed  on 
species  of  atmospheric,  astronomical,  chemical,  and  spectroscopic  interest.  This  work  has  in¬ 
cluded  studies  of  light  asymmetric  rotors  such  as  H2O  and  H^S,  prototype  internal  rotors 
such  as  HOOH  and  CH3OH,  light  asymmetric  rotors  with  electronic  angular  momenta  such 
as  N02,  HOj,  and  HCO,  and  a  number  of  unstable  species  and  ions  such  as  CCH,  NO*,  and 
LiH.  In  addition  to  the  experimental  aspects  of  this  work,  we  have  also  developed  theoretical 
and  computational  methods  for  characterizing  to  microwave  accuracy  these  spectra  over 
wide  spectral  ranges.  Many  (if  not  most)  of  these  studies  were  the  first  to  provide  spectral 
characterization  of  these  species  to  microwave  accuracy  over  large  portions  of  their  rotational 
spectra  in  the  mm/submm.  This  work  is  of  importance  to  the  reported  work  because  the 
basic  spectroscopic  properties  (energy  levels,  transition  frequencies,  transition  moments,  etc.) 
of  the  molecular  species  which  we  studied  must  first  be  well  understood. 

2.  Collisionally  cooled  spectroscopy  at  very  low  temperatures:  As  is  discussed  in  more 
detail  in  Section  HI,  we  have  devoted  a  significant  portion  of  our  effort  recently  toward  the 
development  and  exploitation  of  a  collisional  cooling  technique  which  provides  a  simple 
method  for  the  study  of  gas  phase  processes  at  temperatures  below  S  K.s>*  The  scientific 
motivation  for  this  work  has  been  the  desire  to  study  molecular  interactions  in  the  regime 
where  kT  £  hvr.  To  this  point  we  have  focused  our  attention  on  the  study  of  pressure 
broadening  in  the  region  between  5  K  and  1.5  K  and  on  investigations  of  vibrational  and 
rotational  relaxation  at  77  K.  In  this  regime  we  have  found  that  the  mechanisms  responsible 
for  broadening  make  a  transition  from  being  dominated  by  classical  inelastic  collisions  at 
higher  temperatures  to  a  lower  energy  regime  in  which  resonances  associated  with  quasi¬ 
bound  states  become  important  The  variety  and  richness  of  the  results  that  we  have  obtained 
have  been  beyond  our  expectations. 

3.  Collisional  energy  transfer:  We  have  also  done  considerable  work  on  the  internal 
mechanisms  of  molecular  lasers  and  the  basic  collisional  energy  transfer  processes  that  are 
the  basis  of  them,  especially  optically  pumped  FIR  lasers.7’,  Our  emphasis  has  been  on 
rotational  processes  both  because  this  degree  of  freedom  corresponds  most  closely  with 
mm/submm  devices  and  because  a  number  of  interesting  basic  problems  exist  in  rotational 
energy  transfer. 
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II.  EXPERIMENTAL  AND  THEORETICAL  CONSIDERATIONS 


Because  the  thrust  of  this  report  is  the  study  of  molecular  systems  under  nonambient  con¬ 
ditions,  it  is  useful  to  briefly  consider  the  spectroscopic  consequences  of  significant  changes 
in  temperature  both  on  spectroscopic  parameters  and  the  physics  of  the  systems  being  stud¬ 
ied. 


A.  Oas  Phase  Molecular  Interactions  as  a  Function  of  Temperature 

First  consider  the  temperature  scaling  laws  for  the  interaction  of  gas  phase  molecular 
systems  and  mm/submm  radiation,  For  the  purposes  of  numerical  comparison,  we  will 
assume  that  a  system  originally  at  300  K  is  cooled  to  4  K  and  that  it  may  be  diluted  in  a  colli- 
sional  cooling  gas. 

1.  Line  width;  In  the  mm/submm  spectral  region  experiments  are  typically  run  under 
doppler  broadened  limits.  Since  doppler  broadening  is  proportional  to  Tl/2,  linewidths  at  4  K 
are  typically  an  order  of  magnitude  narrower  and  correspondingly  greater  resolution  and 
measurement  accuracy  are  possible.  In  addition,  since  sensitivity  is  ordinarily  determined  by 
peak,  not  integrated  intensity,  this  factor  under  many  circumstances  also  adds  an  order  of 
magnitude  to  system  sensitivity.  We  have  confirmed  this  expectation  in  a  number  of 
studies.5,4,9'13 

2.  Sensitivity:  Several  temperature  dependent  factors  affect  molecular  absorption 
coefficients.  The  rotational  absorption  coefficient  is  given  by 

a  *  [8rc2v/3ch]  •  l<ml^ln>l2  •  [N/Av]  •  [l-e*hv/kT]  •  [1/Qr] 

where  it  has  been  factored  to  separate  the  individual  temperature  dependent  terms.  The 
behavior  of  N/Av  depends  upon  the  details  of  the  pressure  broadening,  but  simple  theory 
gives 


N/Av  -  T1* 

Qr  is  the  rotational  partition  function  which  approaches  1  in  the  low  temperature  limit  If  we 
consider  molecules  whose  J  »  0  - 1  transitions  fail  at  -  100  GHz,  Qr  is  -  100  - 1000.  The 
term  [l-e'^v^],  which  represents  the  difference  between  induced  absorption  and  emission 
varies  between  unity  in  the  limit  hv  »  kT  and  -1/100  for  transitions  of  v  -  100  GHz  at  300 
K.  Together  these  factors  represent  a  four  or  five  order  of  magnitude  gain  in  absorption  coef¬ 
ficient  for  systems  as  very  low  temperature.  Other  factors  are  more  dependent  on  the  details 
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of  the  system.  The  dilution  ratio  of  the  spectroscopic  gas  in  the  cooling  gas  will  reduce  this 
gain  for  stable  species,  but  the  reduced  pressure  broadening  parameters  associated  with  zero 
dipole  moment  collision  partners  and  the  gains  in  detector  sensitivity  for  low  background  en¬ 
vironments  will  add  to  the  overall  system  gains.  The  latter  factor  will  be  especially 
significant  when  detector  noise  is  the  limiting  factor  in  the  system.  To  give  two  typical 
examples,  the  absorption  coefficients  of  gas  phase  CO  and  HCN  at  4  K  are  25  cm'1  and 
10000  cm*1  respectively.  Thus,  we  conclude  that  even  at  very  large  dilutions  of  the 
spectroscopic  species  in  the  collisional  cooling  gas,  absorption  coefficients  will  be  very  large. 
Again,  this  expectation  has  been  confirmed  in  our  recent  work  and  molecules  with  small 
dilution  ratios  (~  104  •  10  *6)  are  easily  observable  on  an  oscilloscope  screen  in  real  time. 

3.  Physics:  The  spectral  complexity  of  many  systems  is  considerably  reduced  at  low 
temperature.  For  example,  it  has  been  elegantly  demonstrated  by  many  workers  using  free 
jet  expansions  that  the  complexity  of  rotational  structure  in  room  temperature  spectra  can  be 
dramatically  simplified.14'13  In  addition,  rotational  collisional  processes  are  of  both 
fundamental  and  practical  importance  and  many  articles  and  reviews  on  the  subject  have 
been  written.  In  a  review  article,  Brunner  and  Pritchard16  compare  and  contrast  ‘photon’ 
spectroscopy  and  ‘collisional’  spectroscopy.  They  point  out  that  both  spectroscopies  are 
characterized  (in  principle)  by  massive  amounts  of  data  (mostly  redundant)  that  should  be 
calculable  from  a  relatively  small  number  of  ‘spectroscopic*  constants.  However,  the  nature 
of  ‘collisional’  spectroscopy  is  such  that  it  is  much  more  difficult  to  invert  the  problem  so 
that  fundamental  parameters  may  be  recovered  horn  the  experimental  results.  Even  if  these 
parameters  are  given,  it  is  usually  a  formidable  computational  task  to  calculate  the  effects  of 
these  collisions  either  for  the  comparison  of  experiment  and  theory  or  for  use  in  other 
problems.  This  is  because  of  the  participation  of  the  very  large  number  of  thermally 
accessible  rotational  states  in  each  collision  process  and  because  of  the  increased  complexity 
of  molecule-molecule/atom  interactions  in  comparison  with  molecule-field  interactions.  The 
dynamics  of  the  interplay  between  theory  and  experiment  are  further  hindered  by  the  nature 
of  the  available  experimental  data.  Fundamental  parameters  are  recoverable  from  most 
experimental  data  only  via  complex  (and  usually  unrealizable)  deconvolutions.  The  origin  of 
this  difficulty  is  very  similar  to  the  cause  of  the  theoretical  complexity,  the  large  number  of 
thermally  populated  states  and  the  nature  of  the  molecular  interaction.  Thus,  it  is  concluded 
that  the  very  large  reduction  in  the  number  of  thermally  populated  states  at  low  temperature 
will  make  the  relation  between  observables  and  the  fundamental  interactions  substantially 
more  direct.  This  is  one  of  the  principal  motivations  for  our  work.  Finally,  these  very  low 
temperatures  significantly  enhance  the  production  and  lifetime  of  weakly  bound  complexes 
such  as  van  der  Waals  molecules  and  also  provide  additional  opportunities  for  the  study  of 
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ions  and  other  transient  species. 


B.  System  Sensitivity  as  a  Function  of  Temperature 

Over  much  of  the  mm/submm  spectral  region,  the  sensitivity  of  our  experimental 
technique  is  directly  proportional  to  the  sensitivity  of  the  detector  employed.  The  cooled 
detectors  used  in  our  work  approach  being  ideal  detectors  in  that  their  NEP's  approach  being 
limited  by  the  fluctuations  in  the  blackbody  radiation  that  is  incident  on  them  along  with  the 
spectrometer  signal.  In  fact,  if  one  of  these  detectors  is  completely  broadband  and  can  view 
the  entire  300  K  blackbody  spectrum,  its  NEP  is  limited  to  about  5  xlO"11  W/H z1/2.  In  the 
mm/submm  the  use  of  cooled  filters  that  only  pass  the  long  wavelength  portion  of  the  300  K 
blackbody  along  with  the  spectrometer  signal  can  significantly  improve  this,  and  for  a  cutoff 
frequency  of  vm  the  NEP  for  an  ideal  broadband  detector  is  given  in  the  long  wavelength 
limit  by  NEP  ~  kTvm^.17  For  a  vm  •  1000  GHz,  this  correspond  to  -  1045  W.  The  1.5 
K  InSb  detectors  that  we  have  used  in  these  systems  for  many  years2  have  NEP's  in  the 
range  of  1042  *  1043  W/Hz1/2  and  the  0.3  K  3He  detector  mentioned  above  has  a  NEP  of  the 
order  5  x  104S  W/H z1/2,  about  20  times  better.  It  is  our  experience  that  this  has  added  a 
similar  amount  to  our  system  performance.  However,  the  equadon  above  shows  that  the 
construcdon  of  a  more  sensitive  broadband  detector  for  the  mm/submm  spectral  region  is 
not  possible  if  it  must  view  a  room  temperature  experiment. 

There  has  developed  a  substantial  interest  in  the  use  of  the  mm/submm  for  astronomical 
applications.  Because  the  noise  limit  is  this  case  is  often  the  3  K  blackbody  background, 
there  has  been  a  considerable  effort  directed  at  developing  detectors  for  this  low  background 
environment  by  the  use  of  very  low  temperature  techniques.  This  work  has  produced 
detectors  whose  NEP's  are  of  the  order  of  ~1047  W/Hz1/2.  This  is  several  orders  of 
magnitude  better  than  the  limit  set  by  the  thermal  background  for  detectors  viewing  the  long 
wavelength  tail  of  the  blackbody  distribution  and  about  four  orders  of  magnitude  better  than 
the  1.5  K  detectors  commonly  used.  We  are  in  the  initial  stages  of  testing  a  similar  system 
which  is  based  on  an  adiabatic  demagnetization  refrigerator.  Since  much  of  the  proposed 
work  will  take  place  in  an  -4  K  environment,  we  expect  to  be  able  to  realize  much  of  this 
advantage. 
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m.  RESULTS 


A,  Pressure  Broadening  below  5  K  via  Collisional  Cooling 

1.  Experiments]  results:  To  date  we  have  investigated  the  helium  pressure  broadening 
characteristics  as  a  function  of  temperature  for  a  number  of  species  including  CO,  NO,  CH3F, 
HjS,  DC1,  HDO,  OCS,  and  CHjFj  in  the  temperature  region  below  5  K,3'6'9,10'11,12'13  These 
measurements  have  been  made  with  the  collisional  cooling  technique  which  we  have  de¬ 
scribed  in  the  aforementioned  papers.  The  expectation  (based  on  the  calculations  of  sensi¬ 
tivity,  linewidth,  etc.  shown  above)  that  the  collisionally  cooled  technique  should  be  general¬ 
ly  applicable  has  thus  been  borne  out.  The  only  ‘failures’  that  we  have  encountered  have 
been  tractable  to  injector  designs  that  allowed  condensation  or  freezing  for  some  of  the  less 
volatile  of  the  spectroscopic  gases.  Newer  designs  have  completely  eliminated  this  problem 
except  for  HDO,  a  case  where  we  still  have  a  problem  providing  smooth  control  of  the  flow 
rate.  j 

Our  development  of  the  collisional  cooling  technique  was  motivated  by  both  the  general 
ideas  discussed  above  as  well  as  a  desire  to  carry  out  specific  experimental  investigations.  In 
general,  we  wanted  to  be  able  to  experimentally  investigate  a  regime  for  which  hv  £  kT  in  the 
microwave  spectral  region:  thus  having  molecular  systems  in  analogy  to  atoms  at  room  tem¬ 
perature,  with  the  energy  level  spacing  large  compared  to  kT.  However,  only  part  of  the 
analogy  between  molecules  at  very  low  temperature  and  atoms  under  ambient  conditions 
holds  because  at  very  low  temperatures  the  collision  energies  are  small  or  comparable  to  van 
der  Waals  binding  energies.  Thus,  the  molecular  systems  have  the  additional  features  of  sys¬ 
tems  that  can  form  quasibound  states  at  very  low  temperature.1*'20  As  a  consequence,  a  rich 
spectrum  of  results  has  been  observed  that  is  dependent  both  on  the  molecular  species  and  ro¬ 
tational  state. 

Before  showing  experimental  results,  it  is  worthwhile  discussing  some  of  the  qualitative 
features  that  should  be  expected.  Although  it  is  always  dangerous  to  ascribe  too  physical  a 
picture  to  processes  that  are  inherently  nonclassical,  doing  so  provides  in  this  case  a  good 
mechanism  for  classifying  and  understanding  the  results  of  our  experiments. 

In  the  simplest  picture  of  molecular  collisions  and  their  contributions  to  line  broadening, 
the  energy  defect  (which  is  a  measure  of  the  energy  that  is  transferred  between  internal  and 
external  degrees  of  freedom)  and  its  relation  to  kT  plays  an  important  role.  Whenever  this 
defect  becomes  too  large,  the  efficiency  of  the  collisions  in  producing  broadening  and/or  in¬ 
elastic  collisions  is  reduced;  and,  in  the  limit  of  xero  translational  energy,  this  efficiency  also 
becomes  zero.  However,  an  additional  effect  can  become  important  in  very  low  energy  colli¬ 
sions.  Because  of  the  existence  of  shallow  attractive  wells  in  the  intermolecular  potential, 
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resonances  in  both  the  pressure  broadening  and  inelastic  cross  sections  associated  with  quasi¬ 
bound  states  can  exist.  Open  channels  with  significant  strength  result  from  the  contributions 
from  rotational  levels  which  lie  at  low  enough  energy  to  interact  with  the  shallow  wells. 

Thus,  in  the  limit  of  small  collision  energy  there  are  two  countervailing  effects;  the  re¬ 
duction  in  cross  section  due  to  the  growth  in  the  energy  defect  and  the  increase  in  cross  sec¬ 
tion  due  to  the  low  energy  resonances.  Because  both  of  these  effects  are  sensitive  functions 
of  the  intermolecuiar  potential  and  the  energy  level  structure  of  the  molecule,  a  wide  variety 
of  functional  dependencies  of  cross  section  on  temperature  are  observed.  These  qualitative 
expectations  have  been  borne  out  by  diagnostic  calculations  of  Palma  and  Green  for  the  CO  • 
He  system31  and  by  the  experimental  results  described  below, 

These  ideas  lead  to  the  expectation  that  the  pressure  broadening  cross  sections  of  species 
with  the  most  widely  spaced  energy  levels  should  have  the  fewest  resonances  in  their  cross 
sections  and  the  most  rapid  fall  off  in  cross  section  because  of  large  energy  gaps.  As  shown 

in  Fig.  m.A.l,  NO.DCland 
HjS  all  exhibit  monotonically 
decreasing  cross  sections  with 
decreasing  temperature.  How¬ 
ever,  even  in  these  simplest 
cases  the  results  are  far  from 
‘geometric’:  the  cross  section 
at  4.3  K  for  NO  is  modestly 
above  its  300  K  value  of  -22 
A2,  for  DC1  very  near  its  300  K 
value  of  -  10  A2,  but  for  HjS 
much  below  its  300  K  value  of 
-23  A2.  In  a  contrasting  exam¬ 
ple,  shown  in  Fig.  IC.A.2,  CH3F 
(with  its  much  more  closely 
spaced  energy  levels)  has  a 
cross  section  at  4  K  that  is  sub- 


Fifun  m.A.1  Helium  pressure  broadening  croti  sections  of  NO. 
Da.  end  H^S  below  S  K. 
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Figure  m.A.2  Helium  pressure  broadening  cross  sections  for  the  K 
■  0  component  of  the  J  ■  2  •  3  transition  of  CH^F.  The  dotted  lines 
indicate  possible  connections  between  the  classical  region  at  high 
temperature  and  the  cross  sections  below  3  K. 


s tan ti ally  higher  than  its  room 
temperature  value.  The  broken 
line  in  the  figure  is  drawn  to  in¬ 
dicate  the  uncertainties  associat¬ 
ed  with  how  the  300  K  mea¬ 
surement  connects  with  our  re¬ 
sults  below  5  K.  We  have  also 
observed  about  a  25%  variation 
among  different  transitions  at  4 
K  in  CH3P.  It  is  useful  to  note 
that  although  CH3F  (and  other 
species  with  fairly  dense  energy 
level  structures)  show  cross  sec¬ 
tions  at  4  K  that  are  significant¬ 
ly  higher  than  their  300  K  val¬ 
ues,  all  species  observed  to  date 
show  either  flat  or  decreasing 


cross  sections  at  T  -  1 .8  K. 

Figure  1U.A.3  shows  a  comparison  lietween  the  cross  sections  for  the  lowest  rotational 
transitions  of  NO  (J  ■  1/2  -  3/2)  and  CO  (J  *  0  - 1).  As  would  be  expected  for  these  similar 

molecules,  the  cross  sections  at 
room  temperature  are  essential¬ 
ly  identical.12,22  In  addition,  for 
CO  measurements  down  to  77 
K  reveal  an  essentially  un¬ 
changed  cross  section. 

However,  at  4.2  K  the  results 
begin  to  rapidly  diverge,  with 
the  NO  cross  section  rapidly  de¬ 
creasing  and  the  CO  increasing. 

Figure  m.A.4  shows  the 
variation  with  J  of  the  cross  sec¬ 
tions  of  CO  for  the  region 
below  5  K.  Inspection  of  this 
figure  shows  significant  state  to 
state  variation  in  cross  section, 


Figure  m.A.3  A  comparison  of  the  helium  pressure  broadening 
cross  sections  for  the  lowest  rotational  transitions  of  CO  and  NO. 
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Figure  m.A.4  The  dependence  on  J  of  the  crosi  section  of  CO 
below  5  K.  The  smooth  lines  are  representations  of  the  experimen¬ 
ts  data. 


Tsmpsraturs  (K) 


Figure  m.A.5  Helium  pressure  broadening  cross  section  for  the 
\l  ‘  ^  13  transition  of  CH2F2  below  5  K. 


with  the  expected  qualitative 
behavior,  a  faster  decrease  in 
cross  section  for  the  more  wide¬ 
ly  spaced  levels. 

Recently,  we  have  also 
studied  the  3a  313  transition 
of  CH2F2  and  the  results  are 
shown  in  Fig.  m.A.5.  In  addi¬ 
tion  to  this  being  the  first  asym¬ 
metric  rotor  with  a  reasonable 
dense  energy  level  structure  that 
we  have  studied,  it  contains  the 
lower  pump  level  for  the  well 
known  CH2F2  FIR  laser.  Be¬ 
cause  this  level  is  still  signifi¬ 
cantly  populated  at  4  K,  this 
system  is  an  good  candidate  for 
time  resolved  rotational  energy 
transfer  studies  at  very  low  tem¬ 
peratures.  This  will  be  dis¬ 
cussed  in  more  detail  in  Secdon 
IV. 

2.  Discussion:  Numerous 
calculadons  have  shown  that 
cross  sections  are  especially 
sensitive  to  the  intermolecular 
potential  at  low  temperatures.11' 
iwiias  Because  these 
resonances  require  the  existence 
of  at  least  a  shallow  attractive 
well,  the  depth  of  the  assumed 
well  becomes  critical.  Only  for 
the  CO  -  He  system  have  de¬ 


tailed  calculadons  been  pub¬ 
lished.34  The  results  of  these  calculations  for  the  J  ■  0  -  1  and  J  ■  1  -  2  transitions  are  shown 


in  Fig.  m.A.6.  Comparison  of  Figs.  m.A.4  and  m.A.6  shows  basic  agreement  between  the 
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Enr*y,  cm"' 


Figure  m.A.6  Calculated  helium  pressure  broadening  cross  section 
for  the  J  ■  0  -  landJ-  1  *  2  transitions  of  CO. 


experimental  and  theoredcal 
results  around  4  K,  but  at  lower 
temperatures  there  are 
significant  differences.  The 
general  theoredcal  predicdon 
that  the  J  ®  1  -  2  cross  section 
is  smaller  than  the  J  —  0  —  1  is 
borne  out,  but  for  both 
transitions  near  1.8  K  the 
experimental  values  fall  about 
30%  below  the  predictions  of 
theory.  It  is  thought  that  this  is 
a  result  of  the  proclivity  of  the 
methods  used  to  calculate  the 
intennolecular  potential  to  over¬ 
estimate  the  depth  of  the  attrac¬ 
tive  well.23 

For  CO  -  He  the  rapid  rise 
in  the  calculated  cross  sections 
with  reduced  temperature 
results  from  the  thermal  average 
over  a  number  of  strong,  sharp 
resonances  in  the  cross  section 
at  low  energy.  These 
resonances  are  strongly  depen¬ 
dent  on  the  details  of  the 
assumed  potential,  and  Palma 
and  Green21  have  investigated 
this  dependence  in  a  subsequent 
paper.  They  calculated  two 
additional  sets  of  cross  sections 
in  which  the  depth  of  the 
attractive  well  of  the  CO  •  He 


potential  used  in  the  earlier  work34  was  first  reduced  by  one  half  and  finally  completely 


eliminated.  The  results  of  these  calculations  are  shown  in  Fig.  m.A.7.  The  most  important 


result  was  the  demonstration  that  a  reduction  in  the  attractive  well  reduces  the  number,  width, 


13 


*****  CntTft.t m-* 


j  ■ » ■  i  u»  o,co  _ 

^MCnKBUKa.  ^aar-«««— . 

,waon  w°w  jo  “•■  *»  ^££2r  ^ 

cnt  in  those  figures  is  the 


i* 

r 


•^iwi  "  "  •  *  Tt — 

*"^W  (•«•')  H 


Wgrne  HI.A.8  CMlcultted  helium  M 

tr  Poo  OMlttae  ofHDO.  '"'  b,Mtolln«  *»  Wdon,  rOT  Oofttori^ 4,  ^  ,  y 


10'  J01  -0oo< 


structure  in  the  cross 


*«i0„.  that  reeol,  fiom 


Because  of  the  widely 


spaced  levels  and  geometry  of  this  light  hydride  there  is  significant  variation  among  these 
transitions.  In  addition,  these  calculations  show  a  continuous  monotonic  increase  in  cross 
section  above  energies  of  about  20  cm'1.  Of  particular  interest  is  the  fact  that  the  calculations 
show  that  these  resonances  are  spaced  widely  enough  that  they  survive  the  thermal  averaging 
to  produce  significant  structure  in  the  pressure  broadening  cross  sections.  However,  since 
these  calculations  were  done  on  a  1  cm'1  energy  grid,  additional  calculations  will  be  required 
to  show  that  no  significant  resonances  were  missed.  Additional  resonances,  if  present,  could 
smooth  out  the  pressure  broadening  cross  sections  by  filling  gaps. 

The  calculations  on  CO  •  He  and  HDO  •  He  provide  a  basis  for  a  qualitative 
understanding  of  the  experimental  results  for  the  other  systems  for  which  detailed  calcula¬ 
tions  do  not  exist.  At  the  lowest  temperatures,  none  of  the  observed  transitions  show  an 
increasing  cross  section  with  decreasing  temperature.  Indeed,  many  show  a  rapidly 
decreasing  cross  section.  As  would  be  expected  qualitatively,  the  species  with  the  more 
widely  spaced  levels  show  smaller  cross  sections,  and,  within  CO,  the  higher  J  levels  have 
smaller  cross  sections.  These  observations  are  clearly  consistent  with  the  fact  that  at  the 
lower  temperatures  the  size  of  the  energy  defect  increases  relative  to  the  available 
translational  energy  and,  therefore,  insufficient  energy  exists  to  initiate  inelastic  collisions. 
However,  it  is  perhaps  surprising  that  there  exists  only  modest  state  to  state  variation  within 
CH3F  because  the  separation  of  the  states  of  the  several  observed  transitions  from  their 
neighbors  varies  widely.  In  addition,  those  speoies  such  as  CH3F  and  CH^  which  have 
some  level  spaced  £  kT  show  cross  sections  at  4  K  considerably  higher  than  at  300  K,  pre¬ 
sumably  due  to  the  existence  of  the  aforementioned  resonances. 

Finally,  it  should  be  noted  that  the  basic  question  of  the  asymptotic  behavior  of  cross  sec¬ 
tion  in  the  limit  of  zero  temperature  is  open.  Numerical  calculations  for  the  J  •  0  *  1  of  CO  - 
He  indicate  a  very  rapid  increase  at  very  low  temperatures.  Although  the  J  -  0  - 1  experi¬ 
mental  data  show  a  significant  increase  starting  at  4  K,  the  data  at  the  lowest  available  tem¬ 
peratures  appear  to  be  flattening.  On  the  other  hand,  the  calculations  for  J  ■  1  -  2  show  an  ini¬ 
tially  rapidly  rising  cross  section  with  decreasing  energy,  but  a  flattening  at  the  lowest  ener¬ 
gy.  The  theoretical  calculations  for  HDO  were  done  on  a  coarser  energy  grid  and  more 
uncertainty  exists  about  their  low  temperature  behavior. 

Thus,  although  the  general  correlation  between  density  of  states  and  behavior  of  cross 
section  with  collision  energy  may  be  qualitatively  understood,  a  number  of  important  experi¬ 
ments  remain  to  be  done  and  a  special  need  for  additional  calculations  to  strengthen  the  inter¬ 
play  between  experiment  and  theory  exists. 
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B.  Energy  Transfer  in  Polyatomic  Molecules 

Our  work  on  rotational  and  fast  vibrational  energy  transfer  grew  out  of  an  initial  interest 
in  constructing  models  of  the  internal  dynamics  of  FIR  molecular  lasers  which  would  include 
modem  ideas  of  molecular  collision  dynamics.  Although  the  number  of  degrees  of  freedom 
in  these  molecular  systems  might  at  first  glance  seem  discouraging,  we  used  mm/submm 
spectroscopic  diagnostic  techniques  to  establish  experimentally  the  existence  of  pools  of  lev¬ 
els  which  are  in  equilibrium  with  each  other,  thereby  substantially  simplifying  the  resulting 
models.  Because  of  the  very  high  sensitivity  of  mm/submm  spectroscopy  to  nonthermal 
rotational  populations,  this  could  be  done  with  considerable  certainty.  Examples  of  such 
thermal  pools  include  the  equilibria  of  the  vibrational  bending  modes  with  the  R/T  tempera¬ 
ture  in  the  HCN  discharge  FIR  laser,  all  of  the  J  rotational  levels  within  the  unpumped  K  »  1 
state  of  the  12CH3F  OPFIR  laser,  all  of  the  rotational  levels  within  the  unpumped  symmetry 
species  of  the  13CH3F  OPFIR,  etc.24'7'27  This  work  has  provided  a  quantitative  model  of  the 
complex  HCN  discharge  laser  system  as  well  as  significant  advances  in  our  understanding  of 
OPFIR  lasers.  Perhaps  the  most  interesting  practical  result  of  this  latter  work  was  the  dem¬ 
onstration  and  theoretical  explanation  of  the  operation  of  an  OPFIR  laser  in  a  high  pressure, 
small  cavity  regime  that  was  previously  thought  to  be  forbidden  on  rather  fundamental 
grounds.31 

This  work  is  related  to  other  efforts  in  collisional  energy  transfer.  A  number  of  years  ago 
Oka  reviewed  collision  induced  rotational  relaxation29  and  reported  on  his  pioneering  work 
as  well  as  on  the  results  of  others.30*33  That  work  laid  out  basic  ‘selection*  and  ‘propensity’ 
rules  for  rotational  collisional  energy  transfer  and  presented  experimental  verification  for 
some  of  these,  especially  for  those  based  on  spin  statistics.  Since  then,  other  workers  have 
reported  observations  of  infrared-microwave  (ir-mw)  double  resonance  effects.34*36 
Recently,  a  number  of  double  resonance  techniques  have  been  used  for  quantitative  studies  of 
rotational  relaxation  in  polyatomic  molecules.37*42 

From  all  of  this,  it  has  become  clear  that  vibrational  energy  transfer  and  rotational  energy 
transfer  are  not  separate  and  factorable  subjects.  This  is  especially  true  either  when  systems 
have  very  fast  vibrational  processes  which  are  therefore  convolved  with  rotational  processes 
or  when  Coriolis  or  other  mixings  give  rotation  like  character  to  ‘vibrational’  transitions.  A 
nice  example  of  the  latter  is  the  work  on  isotopes  of  formaldehyde  by  Haub  and  Orr.40  In  an¬ 
other  example,  we  have  shown  in  a  quantitative  modeling  of  the  HCN  FIR  discharge  laser 
that  a  similar  Coriolis  interaction  is  the  principal  dissipative  agent  for  the  vibrational  state 
containing  the  upper  laser  level.26 

In  the  work  described  below  time  resolved  mm/submm-ir  double  resonance  methods 
were  used  to  obtain  the  experimental  data  which  formed  the  basis  for  the  development  of  nu- 
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merical  models  which  characterize  the  molecular  energy  transfer.  The  experimental  details 
as  well  os  the  results  of  these  investigations  have  been  published  in  a  number  of  papers.7,1,27, 
43  A  brief  discussion  of  the  diagnostic  techniques  used  for  these  investigations  is  included  in 
Appendix  B. 

1.  Modeling:  As  a  first  goal  in  our  modeling  of  FIR  lasers  and  rotational  and  vibrational 
energy  transfer,  we  have  sought  to  build  the  kinds  of  energy  transfer  models  discussed  above 
that  describe  the  systems  in  terms  of  a  finite  number  of  individual  energy  levels  and  thermal 
pools  along  with  the  energy  transfer  rates  that  connect  them  together.  Although  the  resulting 
models  are  too  complex  for  analytical  solution,  we  have  found  that  they  may  be  implemented 
numerically  on  microcomputers  using  relatively  simple  techniques.  In  this  form  nonlinear 
least  squares  techniques  are  used  to  adjust  the  rates  of  the  model  to  match  the  time  resolved 
data  obtained  from  our  experiments.  At  a  more  fundamental  level  we  have  sought  to  obtain 
rates  that  can  be  related  to  fundamental  physical  processes. 

As  an  example.  Fig.  m.B.l  illustrates  the  relevant  energy  levels  and  pumps  for  l2CH3F 
and  13CHjF.  Even  though  vibrational  relaxation  to  the  ground  state  and  rotational  relaxation 
of  the  ‘hole’  pumped  in  the  lower  vibrational  state  are  important  processes,  they  are  essen¬ 
tially  mirror  images  of  processes  in  the  excited  states  and  for  simplicity  are  not  included  in 
the  figure.  Although  spectroscopically  these  isotopes  are  very  similar,  it  is  important  to  note 
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Figure  m.B.l  A  comparison  of  tha  lymmstry  special,  thermal  pools,  and  Donthamal  J  tttss  within  the 
pumped  X  manifold  for  l2CH3F  and  13CH3F.  Areas  are  reptfsentative  of  the  sixes  of  the  partition  function  of 


that  the  C02  laser  pumps  connect  rather  different  J  levels  and  entirely  different  symmetry 
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species  in  the  two  isotopic  species.  One  measure  of  the  significance  of  these  differences  is 
the  rather  large  difference  in  their  characteristics  as  FIR  lasers.44  As  a  result  the  dynamical 
data  obtained  from  the  time  resolved  mm/submm-ir  double  resonance  experiments  on  each 
species  are  not  strongly  correlated. 

Briefly,  we  have  found  it  possible  to  account  for  a  large  and  diversified  body  of  experi¬ 
mental  data  by: 

(a)  Collecting  all  of  the  rotational  states  of  the  symmetry  species  which  does  not  contain 
the  pumped  transition  into  a  single  thermal  pool,  with  the  relative  population  of  the  levels 
within  the  pool  determined  by  the  Boltzmann  factor  and  with  the  total  population  free  to  vary 
in  the  model. 

(b)  Collecting  all  of  the  rotational  states  of  the  symmetry  species  which  does  contain  the 
pumped  transition  into  a  similar  pool,  with  the  exception  of  the  J  states  in  the  pumped  K 
manifold  which  can  have  an  excess  ‘non thermal’  population. 

(c)  Within  the  pumped  K  manifold,  the  nonthermal  populations  of  all  J  are  treated  as  in¬ 
dependent  variables. 

(d)  Additional  pools  for  population  transfer  to  vibrational  states  are  included  as  required 
by  the  magnitude  of  the  vibrational  excitation.  Typically,  for  the  time  resolved  mm/submm- 
ir  measurements  of  energy  transfer  rates,  only  the  v6  is  required  in  addition  to  the  ground  vi¬ 
brational  state  and  the  directly  pumped  v3.7<>  However,  for  strongly  pumped  cw  laser  sys¬ 
tems  we  have  found  that  population  transfer  to  many  additional  states  profoundly  affects  the 
laser  operation.21 

The  principal  energy  transfer  processes  among  these  pools  and  states  are: 

(a)  The  AJ  *  n,  AK  ■  0  processes;  the  fastest  of  which  is  the  electric  dipole  allowed  AJ 
*  1. 

(b)  The  spin  statistic  allowed  AK  ■  3n  processes,  which  distribute  any  nonequilibrium 
population  according  to  a  thermal  distribution  to  all  other  rotational  states  of  the  same  sym¬ 
metry. 

(c)  A  vibrational  swap  processes  which,  for  example,  effectively  transfer  population  be¬ 
tween  the  A  and  E  symmetry  species  of  the  pumped  excited  vibrational  state. 

(d)  Fast  vibrational  processes  (e.  g.  v3  ->  v6  or  v3  +  v3  ->  2v3). 

(e)  Vibrational  relaxation  resulting  from  wall  collisions. 

2.  Results:  Figures  m.B.2  and  m.B.3  show  examples  of  observed  experimental  results 
from  the  time  resolved  mm/submm-ir  experiment.  Figure  m.B.2  shows  the  very  early  time 
emission/absorption  associated  with  the  nonthermal  population  within  the  K  manifold  which 
contains  the  pumped  level  (for  the  case  of  ^GfyF,  the  K  ■  2).  This  early  time  emission  is 
subsequently  followed  by  rotational  relaxation  which  thermalizes  the  J  states  within  the  K 
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manifold,  causing  the  subsequent  absotpdon  maximum  seen  in  all  of  the  traces.  Finally,  the 
excess  population  in  this  K  manifold  is  relaxed  by  transfer  to  other  K  states  of  the  same  sym¬ 
metry,  via  vibrational  swapping  to  the  other  symmetry  species,  by  vibrational  transfer  to 


Figure  m.B.2  Time  resolved  response  of  rotational  transitions  within  the  pumped  K  manifold  for  l2CH3F. 


other  vibrational  states,  and  ultimately  back  to  the  ground  vibrational  state.  Figure 
m.B.3(upper)  shows  the  arrival  (due  to  the  AK  ■  3n  process)  of  this  population  in  another  K 
state  of  the  same  symmetry,  and  Fig.  m.B.3(lower)  shows  the  arrival  (due  to  the  vibrational 
swapping  mechanism)  of  the  excitation  in  a  state  of  the  other  symmetry  of  13CHjF. 


19 


Although  the  details  of  the 
states  and  pools  are  different, 
the  results  are  largely  indepen¬ 
dent  of  isotopic  species  after  ef¬ 
fects  due  to  statistical  weights, 
etc.  are  removed. 

The  time  resolved  experi¬ 
mental  results  shown  above 
were  fitted  via  a  nonlinear  least 
squares  fit  to  the  parameters 
shown  in  Table  m.B.l.  It  can 
be  seen  that  the  agreement  be¬ 
tween  the  results  obtained  for 
the  two  species  is  very  good  as 

Figure  m.B.3  Time  response  due  to  the  AK  ■  3n  (upper  nee)  and  the  fit  to  the  experimental 
vibrational  swap  Gower  trace)  processes.  data  shown  in  Figs.  IILB.2  and 

m.B.3.  The  different  range  of 

AJ  ■  n  rates  which  was  obtained  in  the  two  experiments  reflects  the  difference  in  the  quantum 


Table  QI.B.  1  Colllaional  Energy  Transfer  Rates  for  CH3F. 
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pected  from  the  similar  energy 
level  structure  and  multipole 
moments  of  the  two  isotopic 
species.  The  AJ  ■  n  rates  show 
similar  trends,  with  the  13CH3F 
being  somewhat  faster  than  the 


^F^F.  The  differences  are  of  the  order  expected  given  that  the  population  starts  in 
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J  •  12,  K  ■  2  and  in  13CH3F  in  J  -  3,  K  *  3. 

It  is  worthwhile  briefly  discussing  which  parts  of  the  observed  time  resolved  responses 
are  most  closely  related  to  the  calculated  parameters.  The  ‘vibrational- swap'  process  tran¬ 
scends  the  spin  statistic  selection  rule  and  effectively  transfers  population  in  the  excited  vibra¬ 
tional  state  from  the  symmetry  species  containing  the  pump  to  the  other  symmetry  species. 

A  very  direct  means  of  observing  this  rate  is  to  simply  measure  the  rate  of  population  arrival 
in  a  pool  state  of  the  opposite  symmetry,  for  13CH3F  K  -  1  (K  ■  3  pumped)  and  for  12CH3F 
K  ■  0  (K  *  2  pumped).  Although  the  details  of  this  process  are  not  well  understood,  Mahan43 
(using  the  results  of  London46  and  Margenau47)  has  shown  that  for  cases  of  near  resonance 
the  cross  sections  for  vibrational  collisional  energy  transfer  can  be  considerably  enhanced  by 
the  presence  of  a  large  vibrational  transition  dipole  moment  such  as  exists  in  CH3F.  The  re¬ 
lated  near  resonance  process  (AE  ■  18  cm*1)  v3  +  v3  -4  2v}  has  a  similarly  fast  rate  in 
CHjF.41  A  mechanism  of  this  general  nature  would  lead  to  a  rate  which  would  be  largely 
isotope  independent.  The  rates  of  the  AK  «  3n  processes  are  also  essentially  the  same  for 
each  of  the  isotopes.  Recently,  we  have  done  an  additional  test  of  the  assumption  that  the  AK 
*  3n  process  is  the  same  (i.  e.  thermal)  for  n  ■  1, 2,  etc.  In  these  experiments  12CH3F  was 
used.  K  ■  2  was  pumped,  and  K  *  1,  K  ■  4,  and  K  ■  5  were  observed.  All  were  found  to 
have  the  same  time  response. 

The  experimental  results  for  the  AJ  ■  n  processes  are  more  complex  and  interesting. 
Gearly,  in  the  limit  of  large  energy  gap  (AE  2  kT)  the  energy  gap  must  play  some  role  in  ad¬ 
dition  to  the  role  played  by  the  size  of  the  change  in  J.  Effects  of  this  type  are  well  known  in 
simple  diatomic  species  although  considerable  controversy  exists  about  'energy  gap’,  'power 
law',  etc.  explanations.49*31  However,  in  both  of  the  examples  here,  the  energy  gap  is  smaller 
than  kT  and  it  might  be  reasonably  expected  that  it  would  induce  relatively  small  differences 
between  the  values  of  the  AJ-n  rates  of  the  two  isotopic  species.  Figure  m.B.4  shows  a  plot 
of  the  AJ«n  rates  obtained  for  each  as  a  function  of  AJ.  The  similarity  of  these  two  sets  of 
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data  seems  to  indicate  that  the 
energy  gap  and  quantum  num¬ 
ber  differences  are  not  too  im¬ 
portant  in  the  regime  of  these 
experiments. 

Because  the  deconvolution 
of  the  time  resolved  responses 
of  these  complex  systems  might 
appear  to  be  rather  complicated 
and  might  provide  results  that 
are  very  model  dependent,  we 
have  also  done  a  number  of 
‘sensitivity*  checks  on  the  nu¬ 
merical  model.  First,  it  should 
be  recognized  that  at  very  early 
times  (i.e.  immediately  after  the 
Q-switched  laser  deposits  the 
population  in  the  upper  pumped  state),  first  order  (i.  e.  one  step)  processes  dominate.  Simple 
considerations  lead  to  the  conclusion  that  the  slope  of  the  absorption  at  t  ■  0  represents  the 
difference  between  the  AJ  «  n  rates  that  populate  the  upper  and  lower  J  levels  for  transitions 
within  K  *  2.  Careful  numerical  investigation  shows,  however,  that  the  definition  of  early 
time  depends  upon  the  size  of  n  because  for  large  n  with  fairly  slow  AJ  ■  n  rates,  faster  two 
step  processes  can  become  important  at  fairly  early  times.  Nevertheless,  we  have  been  able 
to  obtain  from  these  early  time  slopes  and  simple  analytical  expressions  results  that  are  in 
good  agreement  with  the  results  of  the  numerical  models. 

In  addition,  we  have  carried  out  a  number  of  numerical  investigations  of  small  modifica¬ 
tions  of  the  theory  (e.  g.  making  the  AJ  ■  n  rates  functions  of  J;  making  the  AK  ■  3n  rates 
small  functions  of  J,  etc.).  Although  these  can  make  small  changes  in  the  details  of  the  fit  at 
medium  and  long  times,  these  changes  are  modest  and  the  rates  calculated  from  the  fit  are 
largely  independent  of  these  choices.  Finally,  the  good  agreement  between  the  results  ob¬ 
tained  from  independent  analyses  of  the  two  isotopic  species  is  good  evidence  that  the  rates 
obtained  from  the  model  are,  in  fact,  closely  related  to  the  physical  processes  with  which 
they  are  identified. 
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Figure  in.B.4  Plot  of  AJ  ■  n  rates  for  12CH3F  and  13CH3F. 
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C.  Millimeter  and  Submillimeter  Technology 

1.  A  Millimeterwave  Ouasioptical  IMP  ATT  Combiner  Source:  Conventional 
waveguide  power  combiners  are  limited  in  power  output,  efficiency,  and  number  of  sources 
that  may  be  combined  in  the  millimeter  wave  region.  This  limitation  is  a  consequence  of  the 
requirement  that  linear  dimensions  of  conventional  waveguide  resonators  be  of  the  order  of 
one  wavelength  to  achieve  acceptable  mode  separation  and  to  avoid  multimode  operation. 

On  the  other  hand,  quasi-optical  resonators  have  linear  dimensions  large  compared  to 
wavelength  and  they  offer  an  attractive  approach  to  overcome  these  limitations.  In  addition, 
the  Q  of  such  resonators  is  high,  resulting  in  a  pure  output  frequency  spectrum  from  a  quasi- 
optical  millimeter  wave  power  combiner.  An  experimental  investigation  of  quasi-optical 
millimeter  wave  power  combiners  based  upon  the  theory  developed  earlier  by  Mink92  is 
discussed  in  this  paper.  The  approach  utilizes  an  array  of  source  elements  placed  at  the 
surface  of  one  reflecting  surface  on  an  open  resonator.  For  testing  purposes,  energy  is 
extracted  at  the  other  reflector  through  a  conventional  waveguide. 

To  experimentally  investigate  the  feasibility  of  quasi-optical  power  combining  of 
millimeter  wave  sources,  an  approach  which  combines  a  wave  beam  or  Fabry-Ferot  resonator 
is  used  as  the  combining  element  and  IMPATT  (Impact  Avolanch  Transit  Time)  diode 
sources  are  configured  in  a  rectangular  array  with  each  diode  connected  to  a  loop  coupling 
element.  The  cavity  configuration  is  shown  in  the  figure  below.  For  ease  of  fabrication,  a 
wave  beam  resonator  of  circular  symmetry  is  employed  in  this  experiment.  The  resonator 
consists  of  two  reflecting  surfaces  which  are  large  relative  to  the  operating  wavelength.  One 
surface  is  a  planar  reflector  and  the  other  is  curved  such  that  it  tends  to  "focus"  the 
electromagnetic  energy  at  the  planar  reflector. 


Quasi -opticil  cavity. 
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A  goal  of  this  experiment  was  to  show  that  a  large  number  of  sources  could  be  combined 
and  that  the  IMP  ATT  diode  sources  would  cohere  with  each  other  in  the  combiner 
environment.  Also,  for  this  experiment,  it  was  decided  to  utilize  commercially  available 
IMPATT  sources  (Hughes  60  GHz  CW  IMP  ATT).  Because  such  devices  are  already 
mounted  and  since  one  side  of  the  diode  is  connected  to  a  heat  sink,  one  could  not  easily 
employ  the  coupling  technique  of  balanced  electrical  dipoles  as  discussed  by  Mink.32  To 
achieve  efficient  coupling  with  a  dense  array  of  elements  it  was  decided  that  loop  coupling 
could  be  utilized.  Other  techniques  could  have  been  employed  however,  but  they  are  not 
easily  extrapolated  to  a  large  number  of  sources.  Space  limitations  here  prevent  the 
presentation  of  the  theory  of  loop  coupling  which  will  be  published  at  a  later  date.  Bach 
source  is  attached  to  a  a  small  radiating  loop  as  shown  in  the  figure  below.  The  diode  heat 
sink  is  embedded  in  a  copper  rod  1  mm  greater  in  diameter,  A  filter  network  is  then  milled  in 
the  copper  as  well  as  a  slot  to  permit  passage  of  the  loop  wire  to  the  diode  cap.  This 
assemblage  is  then  embedded  in  the  plane  surface  of  the  reflector  and  biasing  currem  are 
passed  through  the  niter  network  and  the  loop  to  the  diode.  The  filter  network  consists  of  a 
series  of  quarter  wavelength  transmission  lines  designed  to  present  a  low  impedance  at  the 
point  where  the  loop  passes  through  the  ground  plane. 


FILTER 


LOOP  SLOT 


Diode  package  and  rod  asaembly,  side  and  front  view. 
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Quasi-optical  resonators  are  based  upon  reiterative  wave  beams  or  beam  modes.  These 
modes  were  first  described  by  Ooubau  and  Schwering53  and  for  the  purposes  of  this 
experiment  only  the  lowest  order  mode  is  of  interest.  In  the  transverse  direction  the  lowest 
order  mode  may  be  expressed  by  a  Gaussian  distribution.  In  directions  transverse  to  the 
direction  of  propagation,  characteristic  dimensions  of  the  fields  within  the  wave  beam 
resonator  are  much  larger  than  those  of  a  conventional  waveguide.  They  range  from  about 
10  to  100  wavelengths. 

As  previously  mentioned,  Hughes  60  GHz  CW IMPATT  diodes  were  used  to  test  the 
theory  of  quasi*optical  power  combining.  Diodes  selected  were  silicon  double  drift  devices 
with  a  doping  profile  of  p*pnn+,  where  the  p+  and  n*  regions  allow  for  ohmic  contacts  to  be 
made  to  the  external  circuit.  As  commercially  packaged,  the  p+  region  is  bonded  to 
metalized  diamond  which  is  thermally  and  electrically  conductive.  The  diamond  is  then 
embedded  in  a  0.3  x  0.4  cm  cylindrical  copper  encasement  to  serve  as  a  heat  sink  and  mount. 
The  n+  region  is  coupled  to  a  gold  band  and  then  to  a  gold  button  0.09  cm  in  diameter  sealed 
to  a  quartz  ring  standoff  0.01  cm  above  the  active  device  as  shown  below.  These  packaging 
dimensions  are  believed  to  be  critical  to  the  experimental  results  discussed  later.  The 
IMPATT  diodes  are  designed  to  operate  at  250°C  with  a  MTBF  of  10,000  hours.  At  the 
rated  bias  current,  approximately  7  watts  of  heat  would  be  dissipated  from  each  diode.  The 
diodes  used  were  certified  individually  by  the  manufacturer  with  the  following  intrinsic 
variations:  capacitance  ~  1.50  to  1.91  pf;  breakdown  voltage  (minimum  voltage  for  current  to 
flow)  -22.4  to  24.6  volts;  bias  voltage  for  rated  CW  frequency  -26.2  to  29.3  volts;  rf  power 
-210  to  280  mW;  rated  frequency  -59.0  to  61.6  GHz.  Differences  in  doping  concentrations, 
device  thickness,  and  operating  temperature  (which  induces  thermal  expansion  in  the  p  and  n 
regions)  affect  the  transit  time  of  the  electrons  and  hence  the  rf  frequency  for  a  particular 
diode.54 
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IMPATT  diode  package  configuration  with 
coupling  loop,  adapted  Bom  Kuno.4 
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The  combiner  design  allowed  for  21  diode  rod  assemblies  spaced  9.2  mm  center  to 
center  in  a  modified  rectangular  pattern.  Each  rod  was  embedded  in  a  4x4x0.5  inch  copper 
plate,  as  shown  below.  The  plate  served  as  the  cavity  plane  reflector  and  an  effective  heat 
sink.  The  diode  spacing  was  taken  from  Mink32,  which  gives  normalized  source  spacing  with 
respect  to  fractional  mode  power  and  driving  point  resistances  for  various  array  sizes.  He 
concluded  that  efficient  coupling  could  be  expected  in  the  transverse  fundamental  mode 
when  the  diodes  were  properly  spaced.  The  ultimate  goal  is  to  combine  many  individual 
diodes  to  obtain  one  high  power  source.  Theory  predicts  that  if  a  single  diode  couples  1  mW 
of  rf  power  into  the  fundamental  mode  of  a  quasi-optical  power  combiner,  then,  with  diode 
spacing  and  cavity  dimensions  held  constant,  a  3x5  array  of  diodes  would  combine  to  yield 
300  mW,  a  7x7  array  -  630  mW,  and  a  9x9  array  -  800  mW.  This  is  satisfying  on  physical 
grounds  because  in  the  limit  of  few  diodes  their  "antenna  patterns"  excite  many  modes  in 
addition  to  the  one  coupled  into  the  output  waveguide.  As  the  diode  array  begins  to  fill  the 
spot  size  of  the  mode  at  the  plane  reflector,  its  "antenna  pattern"  is  then  coupled  primarily 
into  a  single  mode.  Finally,  as  the  array  becomes  larger  than  the  spot  size,  the  power  radiated 
by  the  outer  diodes  is  lost. 
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Planar  reflector  with 
embedded  diode  assemblages. 


The  manufacturer  calibrated  each  diode  in  a  standard  IMPATT  diode  test  circuit  in 
which  a  diode  in  reduced-height  waveguide  is  current  biased  and  rf  exits  to  full-height 
waveguide  through  a  stepped  transition  section.  Such  a  test  circuit  was  obtained  and  a 
sample  diode  was  found  to  operate  at  its  certified  operating  point  discussed  earlier.  These 
sources  operate  in  a  low  Q  mode  and  do  not  exhibit  high  spectral  purity.  Power  is  typically 
distributed  across  a  bandwidth  of  several  hundred  kHz  or  more.  Conversely,  a  quasi-optical 
cavity  will  provide  a  significantly  higher  Q.  Since  the  spectral  purity  of  an  oscillator  is 
related  to  the  Q  of  its  resonant  circuit,  the  possibility  of  significantly  enhanced  spectral  purity 
exists.  This  could  be  dramatically  observed  by  changing  the  cavity  spacing  so  that  the  Q  of 
the  modes  varied.  More  than  an  order  of  magnitude  variations  in  spectral  purity  occurred 
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and,  in  favorable  cases,  HWHM  (half  width  at  half  maximum)  bandwidths  that  approached 
10  kHz  were  observed.  * 

For  a  desired  operating  frequency  of  60  GHz,  the  coupling  loop  area  was  calculated  to 
0.4  mm2.  Teflon  and  enamel  coated  30  to  40  gauge  copper  wire  was  tested.  The  diode  was 
current  biased  over  the  complete  range  up  to  300  milliamps.  This  produced  rf  between  40 
and  45  GHz.  Detectable  oscillations  could  not  be  found  above  45  GHz  regardless  of  bias. 
Bias  current  ranged  from  90  to  150  milliamps  when  rf  was  detected,  half  the  nominal  current 
for  60  GHz  operation.  Output  power  was  disappointingly  low.  The  frequency/current 
relationship  in  the  40  to  45  GHz  range  at  these  low  bias  currents  was  distinctly  linear,  an 
interesting  but  not  a  useful  result  for  this  experiment. 

The  disparity  between  calculated  and  observed  oscillation  frequency  for  the  0.4  mm2 
loop  area  is  theorized  to  be  a  result  of  the  packaging  dimensions  and  the  bias  filter  network. 
The  figure  above  shows  that  beneath  the  gold  cap  is  an  area  in  which  a  gold  band  bridges  the 
length  of  the  quartz  ting  connecting  the  cap  to  the  n+  diode  region.  Furthermore,  the  band 
orientation  beneath  the  cap  is  unknown.  This  unexposed  area  may  contribute  additionally  to 
the  0.4  mm2loop  area  constructed  externally  to  the  cap  and  is  loaded  with  dialectrics. 
Accordingly,  the  loop  slot  was  milled  down  close  to  the  cap  to  reduce  the  total  loop  area. 
Better  impedance  matching  and  increased  radiated  power  were  anticipated. 

The  results  were  encouraging.  It  was  possible  to  achieve  an  order  of  magnitude  greater 
output  power  at  bias  currents  above  200  milliamps  at  a  specific  frequency.  Lost,  however, 
were  the  broad  tunability  characteristics  previously  found  with  larger  loops  and  lower 
currents.  Rather,  the  diode  was  tunable  over  a  700  to  800  MHz  range  for  bias  currents 
between  200  and  270  milliamps.  It  was  difficult,  though  possible,  to  achieve  operating 
frequencies  as  high  as  42  GHz.  This  necessitated  pressing  the  wire  loop  around  the  diode 
cap  to  get  the  smallest  external  loop  area.  Possibly  owing  to  the  variability  of  the  unknown 
area  beneath  the  cap,  most  diodes  more  easily  produced  a  radiation  in  the  33  to  37  GHz 
range.  Refinements  to  the  loop  area  would  produce  output  near  35  GHz  for  most  diodes. 
Number  34  gauge  enamel  coated  wire  yielded  the  best  performance  for  both  output  power 
and  durability.  With  a  Hewlett-Packard  8555  spectrum  analyzer,  high  spectral  purity  was 
observed  even  when  the  cavity  length  was  reduced  to  a  few  inches.  At  the  designed  length  of 
10  inches,  the  loaded  cavity  Q  was  measured  to  be  5000,  For  a  single  oscillating  diode  in 
this  configuration,  spectral  purity  was  measured  at  50  kHz  HWHM  bandwidth. 

Frequency,  bias  current,  and  loop  area  were  defined  for  several  diodes  and  the  quest  to 
phase  couple  two  and  more  diodes  in  the  cavity  was  launched.  Ka  band  (26.5  -  40  GHz) 
waveguide  replaced  the  V  band  waveguide  to  accomodate  the  newly  selected  35  GHz  region 
of  operation.  A  wavemeter,  crystal  detector,  oscilloscope,  and  spectrum  analyzer  were 
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connected  to  measure  the  frequency  and  power  of  the  diode.  Two  diodes  were  selected  and 
their  loop  areas  were  adjusted  for  operations  near  35  GHz.  When  within  1  GHz  of  each 
other,  they  became  candidates  for  coupling  in  the  combiner.  Coupling  with  two  diodes  was 
not  spontaneous,  but  easily  achieved  by  alternately  adjusting  the  bias  current  and  cavity 
length.  Precise  current  tuning  was  achieved  with  precision  potentiometers.  Cavity  lenth  was 
adjusted  with  a  micrometer-driven  translation  table  on  which  the  plane  reflector  was 
mounted.  The  acdvity  was  easily  monitored  from  either  the  dc  level  on  the  scope  or  on  the 
spectrum  analyzer.  Once  the  diodes  were  coupled,  the  spectrum  analyzer  revealed  that  they 
remained  so  despite  varying  bias  current  by  as  much  as  5%.  Coupling  enhanced  spectral 
purity  which  was  measured  to  be  25  kHz  HWHM  bandwidth.  In  combination,  two  diodes 
produced  about  four  times  more  power  than  the  individual  diodes,  in  agreement  with 
theoretical  predictions.  A  third  diode  was  prepared  in  a  similar  fashion  and  bias  tuned  to 
near  the  operating  frequency  of  the  first  two  diodes.  Phase  locking  of  this  third  diode  was 
straight  forward.  It  is  expected  that  the  rest  of  the  array  can  be  filled  in  by  continuing  this 
process. 

In  conclusion,  power  combining  with  millimeter  wave  IMPATT  diodes  in  a  quasi-optical 
cavity  has  been  demonstrated.  Power  exceeding  merely  the  sum  of  the  individual  diode 
outputs  has  been  observed  as  predicted.  The  diodes  oscillating  in  a  high  Q  quasi-optical 
cavity  produce  radiation  of  high  spectral  purity.  Formidable  obstacles  still  exist  in 
understanding  impedance  matching  of  wire  loop  coupled  antennas  to  IMPATT  diodes. 
Inherent  variability  in  commercially  packaged  IMPATT  diodes  has  complicated  the 
theoretical  calculation  of  oscillating  frequencies  due  to  the  introduction  of  unknown  coupling 
loop  geometries.  These  problem  areas  may  eventually  be  overcome  by  epitaxially  growing 
diodes  on  a  single  substrate  and  fabricating  coupling  loops  to  precise  dimensions  with 
microstrip  techniques  to  produce  a  large  array  of  millimeter  wave  diodes  in  a  monolithic 
integrated  circuit. 
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2.  Submillimeterwave  Imaging :  We  have  also  carried  out  at  a  modest  level  a 
submillimeterwave  passive  imaging  project.  The  impetus  for  this  work  has  been  our  interest  in 
the  development  of  very  sensitive  broadband  detector  techniques  for  this  spectral  region.  Since 
these  detectors  are  the  key  to  passive  imaging  in  this  spectral  region,  several  observations  were 
made  both  within  the  laboratory  and  in  the  external  environment.  These  used  small  (4**  -  8") 
optics  and  slow  scan  techniques.  Orders  of  magnitude  enhancements  would  result  from  the  use 
of  larger  antennas  and  focal  plan  arrays.  One  of  the  most  remarkable  results  from  this  work  was 
that  when  we  presented  our  results  at  an  ARO  sponsored  workshop  on  millimeter  and 
submillimeter  wave  imaging,  the  group  in  attendance  concluded  that  they  were  the  only 
submillimeter  wave  images  that  had  ever  been  produced! 

In  the  selection  of  an  active  or  passive  imaging  system,  basic  physical  considerations  impose 
a  strong  wavelength  bias.  At  long  wavelengths  active  systems  (e.  g.  radar)  are  ordinarily 
selected,  while  at  short  wavelengths  passive  systems  are  predominant  (e.  g.  eyesight,  film, 
thermal  sensors,  etc.).  This  is  largely  due  to  the  fact  that  thermal  emission  falls  rapidly  at  long 
wavelength,  peaking  at  lOji  for  objects  of  -300  K  and  at  -6000  A  for  objects  illuminated  by  the 
sun. 

Our  work  was  based  on  theoretical  and  experimental  results  which  showed  that  many  of  the 
useftil  features  of  passive  imaging  at  short  wavelengths  can  also  be  realized  at  much  longer 
wavelengths  (300|i  -  3000|i)  by  the  use  of  appropriate  detector  technology  and  that  systems 
based  on  these  longer  wavelengths  have  a  number  of  unique  and  useftil  properties  including,  but 
not  limited  to: 

(1)  Much  improved  penetration  of  clouds,  fog,  smoke. 

(2)  Different  target  signatures. 

(3)  Much  less  well  developed  countermeasures. 

The  signal  to  noise  ratio  of  a  broadband  thermal  detector  is  given  by 

S/N  m  [PJNEPl  l  At]113 

where  Ps  is  the  power  incident  on  the  detector,  NEP  the  noise  equivalent  power  of  the  detector, 
and  At  the  integration  time.  In  the  mm/submm  spectral  region  it  is  ordinarily  a  good 
approximation  to  assume  that  the  Rayleigh-Jeans  limit  applies  and  that  the  signal  power  from  the 
target  (assumed  here  for  simplicity  to  be  a  black  body)  is  given  by 

P,  -  2*T#  A/ 
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where  k  is  Boltzmann's  constant,  Av  is  the  bandwidth  of  the  detector,  and  Ts  the  temperature  of 
the  target.  The  resulting  signal  to  noise  ratio  is53 


SIN  ■ 


WF 


For  T,  -  300  K,  Av  -  10  GHz,  At  -  10  msec,  and  NEP  -  1043  W/Hz4/2,  a  S/N  -  104  is 
calculated.  This  number  would  be  reduced  for  emissivides  <  1,  temperatures  <  300  K,  but  would 
be  increased  by  bandwidths  >  10  GHz.  In  any  event,  it  is  a  large  number  and  attests  to  the  basic 
sensitivity  of  broadband  detectors  at  long  wavelength.  Thus,  we  conclude  that  state-of-the-art 
thermal  detectors  operating  in  the  mm/submm  can  detect  small  temperature  differences  with 
moderate  time  constants. 

A  laboratory  demonstration  of  this  was  reported36  at  the  12th  International  Conference  on 
Infrared  and  Millimeter  Waves.  The  figure  below  shows  a  representative  image  of  a  man  with 
glasses.  This  was  made  using  a  raster  scanned  system  with  a  30  msec  time  constant,  a  bandwidth 
from  -100  -1000  GHz,  a  detector  NEP  of  -10‘15  W/Hz1/2  and  an  8"  mirror.  The  dynamic  range 
and  signal  to  noise  of  the  image  vastly  exceed  the  grey  scale  capabilities  of  the  display  shown. 
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